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1 Introduction 



The studies [I], [21,13 have discovered a new scaling law for seismic events 
on the phase space location-time-magnitude. Specifically, the distribution 
density for time r between two consecutive events of magnitude m > rric has 
the form 



p^t) = XfiXt) (1) 

where A is the rate of events with m > in a given area G, while / is a 
universal function that is independent of the choice of G and cutoff magnitude 
nic- The relation ([T)) is astonishing, being tested (as it has been in P]) for 
a very wide range of rric (between 2 and 7.5), for seismic regions G of very 
different linear size L (between 20 km and the size of the Earth), as well 
as for different catalogs, both regional and global ones, and different time 
periods. 

The parameterization of / seems not yet to have settled down. According 
to 0: 



f{x)=cx^ ^exp(— x/a) (2) 

in the region x > 0.05 with 7 = 0.74±0.05 and a = 1.23±.15. The pioneering 
work ^ uses a parameterization of (j21) for the whole range of x with 7 = 0.1 
(see refined estimates in 0). This allows the behavior of f{x) about zero to 
be interpreted in terms of the Omori law. 

The subsequent discussion strives to answer the following questions: 

What is the distribution of r in examples of synthetic seismicity? 

What is the physical meaning of the parameters 7 and a in (j21)? 

Assuming the form of / to be universal for r, what should it be? 



2 The Model 

Earthquakes frequently form anomalous clusters in space-time. The largest 
event in a cluster is termed the main event. The events that occurred be- 
fore and after the main event in a cluster are called fore- and aftershocks. 
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respectively. It is assumed in a zero approximation that main events consti- 
tute a time-uniform Poisson process. That assumption is widely employed 
in seismic risk studies. 

Aftershocks dominate clusters both as regards their number and duration. 
Their rate as a function of time is described by the Omori law: 



where to is small. Relation (jS)) holds fairly well during the first few tens of 
days (up to a year) with the parameter p = 0.7 — 1.4 At large times 
the value of p becomes greater, occasionally significantly so, making n{t) 
decay in an exponential manner. Taken on the whole, background seismicity 
and spatial interaction do not allow reliable conclusions to be drawn for the 
Omori law at large times. Cases in which holds during decades are unique 



Following the above description, we consider the following model for seis- 
mic events in time. The spatial and magnitude components of an event are 
disregarded for simplicity of reasoning. Let {xi} be a homogeneous point 
Poisson process on a line with rate A*. It is an analogue of main events. Let 
Nq^cH) be an inhomogeneous point process with rate 6{t) + Xo{t). Here, S 
is the delta function, while the presence of 5{t) means that the event t = 
belongs to Nq. The notation A''o(A) defines the number of events Nq in the 
interval A. We will assume that 



This requirement ensures that the total number of events in No is a.s. bounded. 

Consider the infinite series Nl^"\dt), i = 0,±1,±2, ... of independent 
samples of No- The theoretical process is the sum 



The process A"q that has been shifted by the amount Xi can be associated 
with the cluster of the main event Xi. 



n{t) = ct t> to, 



(3) 
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(4) 



N(dt) =J2No\dt-x^). 
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Our task is to describe the distribution of r between two consecutive 
events in A^. The distribution is uniquely specified, because the process N is 
stationary. It is also easy to see that the rate of N is 



According to 



A = A*(l + A). 



d 



P{T>t) = j^P{N{[a,t])=0} 



(5) 



a=0 



and 



P{N{A) = 0) = exp{-A* j P{No{A - x) > 0) dx}. 

The last relation follows from the fact that the main events are poissonian. 
Since the aftershocks make the bulk of a cluster, we shall assume in what 
follows that Ao(t) = for t > 0. Consequently, 



P{No{A) > 0) = 
Combining the above relations, one gets 



1, if G A 
0, if Ac (-00,0). 



CXJ 

P{r>t) = exp{-\* j P{No{u,t + u) > 0)du- X*t} X 



(6) 



X [PiNoiO, t)=0) + J P{No{du) > 0, Noiu, u + t) = 0)]/(l + A). 



We now describe the behavior of the distribution of r near and oo. 
Statement 1. (a) If cluster duration has a finite mean, fd, then 



P{t >t)= exp (-A*(t + f,0)/(l + A) ■ (1 + o(l)), t ^ oo. 
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(b) Let Ao(t) ~ ct-^-^ t oo where < 6 < 1. Then 

P{T>t) =exp{-X*t-0{f-^))/{l + A), t^oo. (7) 
In other words, one has 

hm In P(r > t)/{Xt) = X*/X 

for a Poisson sequence of main events in a broad class of cluster models. In 
terms of the parameterization of that means that 

a = A/A* = 1 + A. 

With a = 1.23 (as in |3^) the main events make a^^ ^ 81% of the total 
number of events. 

The following regularity conditions should be imposed on A^q iii order to 
be able to describe how the distribution density for r behaves for small t: 

PiNoiu,u + t) > I No{u + t} = 1) = o{l), (8) 

P{No{u,u + t) >0 I No{u} = l,No{u + t}) =o{l), t (9) 

where the notation | denotes conditional probability, and A^ol-s} = 1 means 
that there is an event at the point s. We assume in addition that (jHl), © 
hold uniformly in m > 0. 

That last requirement is no limitation for the case of seismic events, con- 
sidering that the rate of cluster events and time relations between them 
seem to be rapidly decaying over time. The requirements (jH}, Q themselves 
ensure that two very closely lying cluster events are not likely to contain 
another cluster event between them, that is, (jHl), © express the requirement 
of sparseness or repulsion for events that are close in time. It follows from 
the obvious inequality P(A'o(A) > 0) < ENo{A) that © will hold, if 
one demands that 

E{No{u,u + t)\A) = o{l), t^O, 



5 



where A = {No{u + t} = 1) in the case © and A = {No{u} = 1, No{u + t} = 
1) in the case Q. 

Statement 2. If © hold, the probability density for r (provided it 
exists) has the following form as t 0: 

oo 

Prit) = [Ao(t) + j \oiu)\uit) du + A(l + A)]/(l + A) • (1 + o(l)), (10) 



where A„(t) = P{No{t + u, t + u + 6)>0 \ No{u} = l)/6 is the conditional 
rate of Nq after time u given a cluster event has occurred at that time. In 
particular, if \o{t) t oo as t ^ and 

Xuit) < kXo{t), 0<t<e, (11) 

then 

1 < Prit)/\oit) < C, t ^ 0. 

In other words, when (jH)), Q hold, the distribution density for r for small t 
is proportional to the rate of cluster events immediately after the main event. 
The statement is not obvious, since any interevent interval is not necessarily 
started by a main event. 

The proofs of the Statements have been relegated to the Appendix. 



3 Examples 

Examples will now be discussed to be able to understand how far the above 
assumptions are restrictive. 

The trigger model. Historically, this is the first seismicity model to appear 
(see 0). It assumes the cluster process A^o to be poissonian. The model has 
not found acceptance in seismicity statistics, because the likelihood of an 
observed sample in that model is technically difficult to use. This does not 
rule out that the model may be helpful, however. 
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Because increments in Nq are independent, the requirement (jS)), Q has 
the form 



u+t 

P{Noiu,u + t) > 0) = j \o{x)dx = o{l), t^O. 

u 

If Ao(a:;) is a decreasing function, one has 

u+t t 

J \o{x) dx < J \o{v) dv = 0(1) . 

u 

Consequently, the decrease of Xq{x) ensures that (jH)), take place uniformly 
in u. The same property of Ao(x) also ensures (jllj) : 



\u{t) = Ao(M + t) < Ao(t). 

We now are going to refine the asymptotic form of Prit) for small t. 

Let Ao(x) be a smooth decreasing function and Xo{t) = ct~^, t < 1. Then 



p,{t) ~ ct-P + Cit"" + C2, t^ 0, 

where a = 2p — 1 forp > 1/2 and a = forp < 1/2. 
This can be seen as follows. When p > 1/2, one has 



00 1 00 

It = j \o{u)\u{t) = J u~^{u + t)~P du + j Xo{u)\o{u + t)du 

00 1 

00 

= cH^'^P j u-^ {I- uyPdu + const + o{l), t ^ 0. 


When p < 1/2, one has 
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oo 



t — 



I 



Ao(m) du + o{l). 







The self-exciting model. A cluster in this model is generated by the fol- 
lowing cascade process. The first event t = is defined as the event of rank 
0. It generates a Poisson process with rate 7ro(t); its events {t[^^} are as- 
cribed rank 1. The procedure then becomes recursive: each event {tf^} of 
rank r = 1, 2, ... generates a Poisson process of its own which is independent 
of the previous ones and which has the rate 7ro(t — t[^^)- The offspring of a 
rank r event are events of rank r + 1, the events of all ranks constituting the 
desired cluster Nq. 

The process N with clusters as described above is known as the self- 
exciting model [Zj or the epidemic type model The model is rather 
popular in the statistical studies and forecasting of seismicity thanks to the 
fact that the predictable component of has simple structure: 



where the ti are events of N{dt) and At = {N{ds), s < t} is a past of the 
process. 

It is easy to see that the rate A of the process N is bounded, if 



Statement 3. (a) The cluster rate function for the self-exciting model is 



E{N{t + 6)- N{t) >0\At) =J2''o{t-ti)-6 + X*6, 



oo 








also, 



A = A^/(l - A^) and A = A7(l - A^). 



AoW = 7ro(t) + TTo * 7ro(t) + TTo * TTo * 7ro(t) + . . . , t>0, 



(12) 
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where * denotes the convolution. 

Let 7ro(t) be monotone near 0, where 7TQ{t) ~ At'"^, < p < 1. Then 

AoW/vro(t) ~ 1, t^O. 
Let vro(t) be monotone at oo, where vro(t) ~ Bt^^^^, < 9 < 1. Then 

(b) The distribution density for t as t ^ has the form 



Prit) = [(1 - A^)Ao(t) + j Ao(x)Ao(x + t)dx + X]-{l + o(l)), t 0. (13) 



Let 7ro(t) be monotone near zero, where 7r(t) ~ At^'^, < p < 1; let 7[Q{t) < 

oo 

where (p is a smooth function, j ip{t) dt < 1, (p{t) ~ ct^^^^, t oo, 



< < 1. Then 

Prit) = 0(Ao(t)) as t i 0. 

r/ie time-magnitude self-exciting model. The self-exciting model is fre- 
quently considered on the time- magnitude space as follows (see, e.g., jH]): 
each event tj (both when a main or a cluster one) is ascribed a random 
magnitude m^. The mj are independent for different ti and have identcal 
distributions with density p{m). The generation of clusters is that described 
above, the only difference being that an event (s,m) generates a cluster 
with rate g(m)7r(t — s). It can be assumed without loss of generality that 
J q{m)p{m) dm = 1. This normalization preserves statements 1, 3 for the 
self-exciting process (t, m) as well, independent of the choice of p{m) and 
q{m). The function Ao(t) as given by (fT^ then corresponds to the cluster 
rate when averaged over magnitude m. For purposes of seismology, p{m) cor- 
responds to the normalized Gutenberg- Richter law, p{m) = pQ-Pi"^-"^o) ^ ^ ^ 
mo while q{m) = e"("^~™'o)(l — a/P) is proportional to the size of the cluster 
that has been triggered by an event of magnitude m. 
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4 The unified scaling law 



According to 3J, the distribution of r depends on the single parameter A, 
see (0). The parameter A is specified by the choice of the area and cutoff 
magnitude rric. This allows variation of A in a very wide range. Experiments 
which test in 3j concern both the Earth as a whole and smaller or larger 
areas of it. One can always select such areas in which seismicity is weakly 
interdependent. For the theoretical analysis of the unified scaling law (0) one 
may be interested in the following 

Statement 4. Assume that it is possible to choose two regions Gi and G2 
with independent stationary sequences of events Ni{dt). If the unified scaling 
law (H)) holds for Gi, G2 and Gi U G2 and f{x) < cx~^ , < 6' < 1 for small 
X, then f{x) = exp(— x). 

Proof By ®, 

Pr{t) = Jg^P{N{0,t)=0}, (14) 

where A is the rate of N{dt) in the region. In virtue of (0) 

pAt) = A/(At). 

Equation ()14|) and the initial conditions for P{N{0, t) = 0} = u{t) having the 
form m(0) = 1 and 'u'(O) = —A specify u(t) uniquely and yield u{t) = (p{\t), 
where 

t 

^{t) = l-t + j{t~ s)f{s) ds (15) 


Since Ni{dt) and N2{dt) are independent, 

P{N{0,t) = 0) = P(iVi(0,t) = 0)P{N2{0,t) = 0), 

where N = Ni + N2 is the sequence of events for Gi U G2- It follows that for 
any t > one has 
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or 



¥^((Ai + A2)t) = v^(Ait)(^(A2t). 



m = Hpt) + Hqt) (16) 

where ip(t) = \rnp(t), p = Ai/(Ai + A2) andp + g = 1. Iteration of (16) yields 
for p = g = 1/2 



or 



^(t) _ jjjent) , . 

m Hen)' ^ ^ 

By ip{0) = 1, ^'(0) = -1 we have ^{0) = and ^'(0) = -1. Using 
L'Hospital's rule we will have 



hm = hm ^^^^"^^^ = t. 

By (fT7|) one has ifjif) = —at or ip(t) = exp(— at). However, in that case 
f{x) = o?e~°'^ and J /(x) (ix = 1, whence a = 1. Statement 4 is proven for 
p = g = 1/2. 

In the general case, p 7^ 1/2, the iteration of (fT^ yields 



k=0 

where ek,n = p^q^''^. 
As above one has 
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with Skn = 0(1) as n — > 00. 

Using (fT^ and the apriori bound 

f{x) < cx~^, < 6 <1, < X < e 
it is easy to show that 

\6,,n\<kf[max{p,q)r^'-'\ 

Therefore we have again ilj{t) = at because 6k^n are small uniformly in k. 
The proof of the statement 4 is complete. 

5 Conclusion 

We have presented a theoretical analysis of the distribution of interevent 
interval r in a point process. It is shown that, when assumptions natural to 
seismic events have been made, the distribution of r may be a function of a 
single parameter, the rate A, provided the distribution is exponential. This 
contradiction means that the nature of the empirical unified distribution of 
r is more complicated. One has to sepatate universal properties of r from a 
visual artifical effects. 

Making very general assumptions, we have found how the distribution of 
T behaves near and 00. As was to be expected IJ, the behavior is related to 
that of the Omori law near zero and to the Poisson character of main seismic 
events, when one deals with asymptotics at infinity. It is these asymptotics 
which essentially make the probability density for r "universal" in (Ij, when 
plotted on a log-log scale. 

The parameterization of the distribution of r put forward in ^ for Xt > 
0.05 has the form f{x) = cx'^~^ exp{—x/a). It was shown above that the 
parameter 1/a can be treated as the fraction of main events among all seismic 
events. The estimate a = 1.23 derived in yields ~ 80%, which can 
hardly be a universal constant. The main events in Italy are 60% among the 
m > 3.5 events (see pUj). 

The factor x"''^ is missing in the formula for f in the models consid- 
ered above. This factor may be replaced (see ((Tj)) by a factor of the type 
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exp(— ca;^~^), if the aftershock rate decays as a power function t~^~^ , < 9 < 
1; the factor degenerates to a constant for > 1. Consequently, it remains 
an open question as to what is the physical meaning of 7. 
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Appendix 



Proof of Statement 1. 

We are going to find the asymptotics of P(r > t) as t — * oo using ((Tj). To 
do tliis, tfie following three limits should be found as t ^ oo: 



A ■= P{N{0, t) = 0) ^ P(iVo(0, oo) = 0) = P{t,i = 0). 
Here, Td is the cluster duration in Nq. 

oo oo 

B:= j P{No{du) > 0, Noiu,u + t) = 0) ^ j P{t,i E du) = P{r,i > 0). 



Consequently, A + B ^ 1. It remains to find the limit for the expression 
under the exp sign in ((7j). One has 

oo oo 

C : = j P{No{u,t + u) > 0)du-^ j P{No{u,oo) > 0)du = 





oo TcJ 

E I lNoiu,oo)>o du = E / du = ETci. 



We have used the notation 1^ : 1^ = 1, if A is true and 1^ = otherwise. 
It remains to substitute the resulting limits into 



P(r >t)= exp{-A*(t + C)}[A + B]/{1 + A). (18) 

We now are going to prove the second part of Statement 1. Let Erd = oo. 
The asymptotics of C then calls for refinement. One has 



t+u 

P{No{u,t + u) > 0) < ENo{u,t + u) = / Xo{v)dv. 
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If \q(v) = c9{l - d)v-^^+^^ for > 1, then 



oo t+u 

C < j du j \Q{v)dv = ct^-\l + o{l)), t^oo, 

u 

as follows from L'Hospital's rule. Relation (|7j) stands proven. 
Proof of Statement 2. 

The distribution density for r can be found by differentiating (jl8|) . To 
do this, we make differences for the functions A, B, C in (fTHj) . One has for 
small S: 

[A{t) - A{t + 6)]/6 = P{N{t,t + 6)>0)6-\l~P{No{0,t)>(}\No{t} = l) 
^ Ao(t)(l + o(l)). 

The last conclusion follows from the regularity requirement imposed on (jSJ. 
Similarly, one has for B{t) using (0): 



[B{t) - B{t + 5)]/5 = j [P{No{du) > 0, No{u + t,u + t + 5) > 0) - 



-P{No{du) > 0, No{u, u + t) > 0, No{u + t,u + t + 6) > 0)] = 

oo 

= J \o{u)K{t) du{l + o{l)), 



where A„(t) = P{No{u + t + 5) - No{u + t) > | No{u} = l)/6. One has 
C'{t) = -A(l + o(l)) for C{t) using ® . It remains to differentiate ^ and 
then to substitute the resulting asymptotic expressions for the derivatives A', 
B', a and the values A(0) = 1, 5(0) = A and C(0) = 0. 
Proof of Statement 3. 

It follows from the description of the cascade generation of A^^o that its 
rate Ao(t) satisfies the integral equation 



Ao(t) = y 7ro(x) Ao(t - x)dx ^ 7ro(t), 





(19) 
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where 7ro(t) is the rate of rank 1 events. Iteration of (|19|) then yields 

Ao(^) = 7ro(t) + TTo * 7ro(t) + TTo * TTo * 7ro(t) + . . . 

If one passes to the Laplace transform, A ^ A, then both relations for Xo{t) 
are reduced to the form 

Ao(.)=^?Gs)/(l-7r(.)). 

Let 7r{t) be monotone near and oo. Assume also that 7r(t) behaves 
like a power law: 7ro(t) ~ Cot~^, t -C 1 or 7ro(t) ~ Cit~^~^, t 3> 1, where 
< p, 6' < 1. In that case the use of the Tauberian theorems (see [H], Ch. 
13 and Ch. 17, §12) yields conclusions of the form Ao(t)/vr(t) — > const as 
t ^ or t ^ oo, respectively. 

We now are going to prove (fT^ . Consider the rate of a pair of events 
in an A^'o cluster: Xo{u,v) = P{N(){du) = l,No{dv) = l)/{dudv), u < v. 
Recalling that this is a cascade generation of Nq, the states u and v in No 
can be derived in two ways. The one is when u and v have no common 
parent except t = 0; the second is when u and v have a common parent z 
in the first generation (a state of rank 1). If the common parent z for u 
and V has rank r > 1, then the probability of that event will be of order 
0{{dz)'^ dudv), which is negligibly small compared with 0{dzdudv). This 
consideration leads to the following equation for Ao(m, v): 

u 

Ao(u, f ) = Ao(u)Ao(f ) + j TCo{z)\o{u — z,v — z) dz, u<v. (20) 



Put at{u) = Xo{u, u + t), bt{u) = Xo{u)Xq{u + t). Then gives 

at{u) = ht{u) + TTo * at{u). 

Whence 

at{u) = ht{u) + ht{u) * (tTo + TTo * TTo + TTo * TTo * TTo + . . .) = 

= bt{u) + bt{u)*Xo{u). (21) 
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We are interested in the conditional rate in a No cluster: 



A«(t) = \o{u,u + t)/Xo{u). 
One has from (PT|) using the notation at and bt: 



Xu{t) = Xo{u + ^) + y ^o{x)Xo{x + t)Xo{u - x) dx/ Xo{u). 



It remains to substitute that expression in ()10|) . One has 



Prit) = [Xo{t) + J Ao(m)Ao(m + t)du-{l+A) + 



+ A(l + A)]/(l + A)-(l + o(l)), t^O. 

However, (1 + A)~-^ = 1 — X-,,, so (fT^ is proved. 

In order to have Prit) = 0{Xo(t)) as t ^ , one has to show that 
Xo{u + t) < kXo{t) for small t. To do this, we demand 7ro(t) = Cit~^, < t < e 

oo 

and 7ro(t) < fit). Here, is a smooth function, f (p{t) < 1 and ip(t) ~ ct^^~^, 



t ^ 1 with < 6 < 1; also, if = tcq for t < e. Then 

is a smooth function. One has Ao(t) < X^{t) in virtue of (fT^ . since ttq < (/J. 
One has X^{t)/(p(t) ^ c as t — *• oo from the power law behavior of at oo 
(see Statement 3a). One also has Ao(t)/7ro(t) ^ 1 as t ^ 0, hence Xo(t) | oo 
ast 0. Consequently, max Ao(t) will coincide with A<^(to) = Ao(to)(l+c>(l)) 

t>to 

starting from some small to- Hence 



Xo{u + t)< max Ao(f ) < max X^{v) ~ Ao(to)(l + o{l)). 

V>to V>to 
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